Finite-temperature conductance signatures of quantum criticality in double quantum 

dots 
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We study the linear conductance through a double-quantum-dot system consisting of an interact- 
ing dot in its Kondo regime and an effectively noninteracting dot, connected in parallel to metallic 
leads. Signatures in the zero-bias conductance at temperatures T > mark a pair of quantum 
(T = 0) phase transitions between a Kondo-screened many-body ground state and non-Kondo 
ground states. Notably, the conductance features become more prominent with increasing T, which 
enhances the experimental prospects for accessing the quantum- critical region through tuning of 
gate voltages in a single device. 

PACS numbers: 73.63.Kv, 73.43.Nq, 72.15.Qm 



Quantum phase transitions (QPTs) occur in the 
zero temperature (T = 0) phase diagram of a sys- 
tem at points of nonanalyticity of the ground-state 
energyiii^ QPTs underhe many fascinating phenomena 
in strongly interacting condensed matter, including the 
metal-insulator transition in disordered systems^ the de- 
struction of antiferromagnetism with doping in high- 
temperature superconductor parent compounds,^ the 
magnetic-field-driven superconducting-insulator transi- 
tion in disordered superconductors,^ and quantum Hall 
plateau transitions.^ Study of most of these QPTs is hin- 
dered by the need to fabricate controlled series of samples 
at different stoichiometrics and/or disorder levels. 

By contrast, it is increasingly apparent that systems 
of quantum dots offer possibilities for exploring QPTs 
(strictly, boundary QPTs involving only a subset of the 
system degrees of freedom) within a single sample. Ad- 
vances in system fabrication, precise characterization, 
and the near suppression of dissipative and incoherent 
environments^, have enabled beautiful experiments on 
multi-dot devices. * This leap forward in experimental ca- 
pability has also spurred much theoretical activity, in- 
cluding several predictions of QPTs in quantum dots in 
the Kondo regime i^ The feasibility of realizing nontriv- 
ial many-body states has been confirmed by the recent 
experimental demonstrations of a two-channel Kondo 
regimei^ and of a singlet-triplet QPTjii 

This Letter predicts robust signatures of QPTs in 
the finite-temperature conductance through a double- 
quantum-dot (DQD) system. A smaller dot ("dot 1") 
exhibits Kondo physics, while a larger dot ("dot 2") is 
effectively noninteracting and lies near a transmission 
resonance. When the dots are connected in parallel to 
external leads, and the system is fine-tuned via applied 
voltages that determine tunneling barriers and the en- 
ergies of individual dot orbitals, a pseudogap in the low- 
energy effective hybridization between dot 1 and the leads 
gives rise to a pair of continuous QPTs between Kondo- 



screened and non-Kondo ground states.— We describe 
how the system can be steered into the vicinity of a QPT 
by monitoring the linear conductance while changing just 
two gate voltages. 

Experimental detection of QPTs necessarily relies on 
finite-temperature manifestations of the underlying T = 
transition. We show that the signatures of quantum 
criticality in the present DQD system become more pro- 
nounced as the temperature is increased from absolute 
zero, a trend that contrasts with the typical behavior 
near an impurity QPT,^ Their temperature dependence 
also allows these signatures to be distinguished from 
other conductance features in the same system. 

Model and conductance calculation. — Consider a DQD 
device in which dot 1 is in an odd-electron-number 
Coulomb blockade valley, and dot 2 has a single level 
near the Fermi level and is effectively noninteracting^^ 
The dots are coupled to left (L) and right (i?) metal- 
lic leads and to each other via tunneling barriers. This 
device is described by a two-impurity Anderson Hamil- 
tonian: 



i? = ^ e^ni^ + Uin^uii + ^ ( A aj^aj^ + H 



Y^ £ik cL^c^ka + Yl (^'^ "-L^eka + H-C.) , (1) 



i,e,k., 



where aj^ creates a spin-cr electron in dot i {— 1,2), n^o- = 

'Ha'^ia' ^^^ 'ika crcatcs a spin-cr electron of wave vector 
k and energy e^k in lead i (= L,R). We assume for 
simplicity that each lead has a density of states p(w) = 
poQiD— |a;|), symmetric about the Fermi energy {lu — 0), 
and that dot-lead couplings are local. We further assume 
that all couplings are real and the device is tuned to left- 
right symmetry, so that we can write Vie — Vi/\/2. 
The linear conductance at temperature T for this DQD 



setup can be obtained from the Landauer formula as: 

/oo 
<h:{~df/dLo)[~luvT{oj)], (2) 
-OO 

r (a;) = 2^po E ^^*L Gy (^) V,r, (3) 

where 50 = 2eV/i, /(w/T) = [exp(u;/T) + 1]-^ is the 
Fermi-Dirac function, and aU Gy(w) in Eq. ([3]) are 
dressed Green's functions, fully taking into account the 
electron-electron interactions on dot 1. 

The standard equations of motion uj{{A:B))^ — 
{{A,B}) = {{[A,H]:B))^ = - {{A; [B , H])) ^ for the re- 
tarded Green's function {{A;B))^ = -i/^ die*"*({A(i), 

5(0)}) allow one to re-express Gij{uo) — {{a^^;aj^))i^ in 

terms of Gu and the bare Green's function G22 , which 
describes the noninteracting dot 2 in the absence of dot 
1. In the wide-band limit |w| <C D^ Eq. ([3]) becomes 

T{uj) = AiGii(w) + 2Ai2[G^2^(c^) (A - ^Al2) Gii(a>)] 
+ A2[l + G^^^H^) (A - ^Au)' Gn{io)] G^^^^c^), (4) 

where Ai = npoV^, A12 — 7TpoViV2, and G22 (w) = 
(lu - £2 +i/^2)~ ■ 

The dot-1 local Green's function Gii{uj) entering Eq. 
([4]) can be obtained-^ by mapping the Hamiltonian ([T]) to 
an effective model of a single dot connected to the leads 
via a nonconstant hybridization function 



A(w) = 7Tp2{uj) [X+iuj- £2) VA1/A2J , (5) 

withp2(w) = A2/{7r[(tj-£2)^ + A|]}. We solve this effec- 
tive model using the numerical renormalization group. ^* 
At T > 0, we compute the spectral function Aii{u!) — 
— 7r~^ImGii(aj), and hence obtain G[i{uj) = ReGii(w) 
via a Kramers-Kronig transformation. At T = 0, where 
Eq. (m involves only Gii(O), it is possible to calcu- 
late G[i{0) directly. All results shown are for Ui — 
0.5D and A2 = 0.02D with temperatures in units of 
Tko = 7.0 X lO^^D, the Kondo temperature in the ref- 
erence case where dot 2 is decoupled (A = A2 = 0) and 
Ui = -2£i = 0.5D, Ai = 0.05D. 

To facilitate interpretation of the results, we note that 
— ImT(w) entering Eq. ^ can be expressed as 

-ImT(cj) = [l-27rA2/92(w)]7rA(w)Aii(w) + 7rA2P2(w) 
+ 27r(cj-e2)A(cc;)p2(w)G;i(cj). (6) 

The term 7rA2P2(^) represents bare transmission through 
dot 2 in the absence of dot 1, and for T ^ A2 yields a 
conductance contribution g2 — 5oA2/(£2 + ^2)- ^^ most 
cases of interest, the term involving G'n turns out to be 
negligible. If, as we assume, the dot-1 level is off reso- 
nance (i.e., |£i| ^ T, Ai), then dot 1 appreciably influ- 
ences g only in the Kondo regime T < Tr- where Aii{uj) 
exhibits a many-body resonance at the Fermi level; the 
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FIG. 1: (color online) (a) Conductance g vs £2 at six tempera- 
tures for a parallel DQD device [inset in (b)] with ei = — C/i/2 
and Ai = 0.05-D. (b) Transmission through dot 1 sets in 
below the Kondo temperature Tk (defined as in[l3). which 
vanishes as £2 — > 0. 



sign of the resulting conductance term gi depends on 
that of 1 — 27rA2P2(^) in the range \ijj\ < 0{T) that de- 
termines g{T). For |£2| 3> A2 (|£2| ^ A2), gi is positive 
(negative) at low temperatures, leading to constructive 
(destructive) interference with 172 ■ 

Tuning to the pseudogap regime. — When the level en- 
ergy in dot 2 is set to £2 = A-\/A2/Ai, the dot-1 
effective hybridization [Eq. ([5|)] vanishes at the Fermi 
level as A(a;) oc w^. The pseudogap Anderson impu- 
rity model, in which A{uj) ex \uj\^ for |a;| —^ 0, exhibits 
Kondo and non-Kondo ground states separated by QPTs. 
Whereas previous theoretical worki^ has focused on ex- 
ponents < r < 1, the proposed DQD setup offers 
a controlled realization of the case r = 2, which fea- 
tures a pair of QPTs. For simplicity, we focus in the 
remainder of the paper on configurations in which the 
dots are connected to the leads purely in parallel, i.e., 
A = [see inset in Fig. [Ijb)]. Then Eq. ([5]) reduces to 



A(w) = Ai(w - £2)V[(^ - £2)' + A 

In order to probe the QPTs, the pseudogap in A(w) 
must be centered on the Fermi energy. Operationally, 
this can be accomplished by tuning £2 (via a plunger 
gate voltage on dot 2) to reach a maximum of g. Figure 
[Ija) illustrates g vs £2 at six temperatures, for fixed £1 = 
-Ui/2 and Ai = 0.051?. For A = 0, A(0) = when dot 
2 is exactly in resonance with the leads: £2 = 0. The 
choice of £1 = — f/i/2 makes £2 = a point of particle- 
hole (p-h) symmetry, and ensures that for |£2| ^ A2 and 
T < Tk, 7rA(0)Aii(0) ~ 1;^^ then, since 7rA2P2(0) ~ 
1, gi almost completely cancels g2- Figure [TJb) shows 
that the temperature range < T < Tk{s2) of the low- 
conductance regime shrinks rapidly as £2 -^ 0. For the 
special case £2 = 0, the pseudogap in A{lu) prevents the 
formation of a Kondo state (effectively, Tk = 0), and 
transport takes place solely through dot 2. At T = 0, 
the resulting conductance exhibits a discrete jump from 
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FIG. 2: (color online) Behavior near the e^^ transition in a 
parallel DQD device with Ai = 0.05D and £2 = 0: (a) Curves 
1-6 show the dot-1 spectral function Aii{uj) for the values 
of Aei = El — ef^ indicated by the corresponding arrows in 
(b). The frequency lo* of the quasiparticle peak in Aii{u}) 
is proportional to Aei. (b) Conductance g vs Aei at six 
temperatures, (c) g vs T/uj* at three values of Aei on the 
Kondo side of the transition. 



(; = for |£2| — * to 5 = go for £2 = [dashed line in 
Fig. [TJa)]. However, this spike broadens at T > into a 
smooth peak rising to g{e2 = 0) ~ go- 

For a general ei ^ —Ui/2, transmission through dot 1 
is still blocked when A(0) = 0. This leads to an asym- 
metric peak in g{e2) at g(0) ~ go — a feature that again 
broadens with increasing T/^ offering a practical method 
for tuning the pseudogap to the Fermi level. 

Tuning to a QPT. — With £2 held at zero, the level 
energy £1 can be varied via a plunger gate voltage on dot 
1. A pair of QPTs, related by p-h duality and located 
at £1 



'-Ic 



'Ui/2 ± |A£ic|, bound a local-moment 
(LM) regime Si^ < £1 < £]^ in which the net spin on 
dot 1 is unscreened at T = 0. Close to either QPT, 
Aii{u!) contains a quasiparticle peak centered ai co — lu*, 
where w* oc £1 — £J^ [Fig. [UJa)]. The peak sets in below a 
crossover temperature ~ |w*|, which on the Kondo side 
is proportional to Tk- This feature in ^11(0;) leads, via 
Eqs. ^ and (O, to a conductance contribution gi < 
that is greatest in magnitude when \uj*\ ~ 4T. Since the 
dot-2 contribution 32 — go is independent of £1, g vs £1 
isotherms [e.g., see Fig.[2Jb)] show a dip at |£i — s^J ex T 
on cither side of a maximum at £1 = Si^. 

It is striking that at T = 0, the conductance shows no 
feature as dot 1 passes through a QPT. At finite tem- 
peratures, by contrast, the DQD device can be tuned to 
the transition by seeking a local maximum in g vs £1. 
This maximum has the identifying characteristics [Fig. 
[DJb)] that the minima on either side are equidistant in 
£1 from £^^, but the dip in g is roughly twice as deep 
on the Kondo side, reflecting the greater weight of the 
quasiparticle peak in that regime. For the parameters 
shown in Fig.[5l^b), the conductance peak becomes more 
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FIG. 3: (color online) (a) Conductance g vs ei for a side-dot 
device (inset) with £2 ~ 0, both for T = at various A values 
and for A = 0.03-D at the labeled temperatures, (b) Phase 
shifts rjii (filled circles) and 7722 (open circles) for A = 0.03D, 
T = 0; g vanishes when sin 7722 = [Eq. ^]. 



prominent with increasing temperature up to T ~ 3Txo, 
and a peak in g remains discernible up to the relatively 
high scale r~6TKo. 

The form of g vs T at fixed £1 is more complicated since 
gi and 32 can have temperature variations of comparable 
magnitude. Figure [2)Jc) shows that in the Kondo regime, 
the peak in |(;i(T)| contributes a shoulder around T = 
\ll!*\ to the overall downward trend dictated by g2{T). 
Similar behavior holds in the LM regime (not shown). 

Differentiating QPT and other conductance features. — 
Conductance peaks similar to those shown in Fig. EJb) 
can also arise, not from proximity to a QPT, but rather 
from interference between a conventional (metallic or 
r — 0) many-body Kondo resonance on dot 1 and a 
noninteracting resonance on dot 2. In experiments, the 
mapping between the gate voltages in a real device and 
parameters of the effective Anderson model will not be 
known a priori. It is therefore important to be able to 
identify unique signatures of a QPT in this system. We 
show below that the temperature dependence of the con- 
ductance peaks serves this purpose. 

For simplicity, we consider the "side-dot" regimei^ 
Ai = in which dot 1 is connected to the leads only 
via the noninteracting dot 2. In this geometry [inset 
in Fig. Ela)], the effective dot-1 hybridization function 
A(a;) ~ ■n}?p2{iS) [from Eq. ([5])] is a Lorentzian of width 
A2 centered at w = £2. Since A(cl') has no pseudogap, 
there is no QPT. 

Figure [3] plots the variation of the conductance with 
the position £1 of the energy level in the side dot 1 while 
dot 2 is held in resonance, i.e., £2=0. At T = 0, 
the conductance drops to zero as £1 approaches the p-h- 
symmetric point £1 = — 1//2, independent of the dot-dot 
coupling A. This can be understood by noting that for 



Ai = and T = 0, Eqs. ^ and ^ reduce to 

<?(r = 0) = -5oA2|G22(0)|sin7722, (7) 

where rju = argGu^O) is the Fermi-energy phase shift of 
electrons scattering from dot i. Figure [3I^b) shows that 
in a window about the p-h-symmetric point, the dot-1 
phase shift exhibits a plateau rjn ~ — 7r/2 characteris- 
tic of the Kondo state. ^^ This additional phase shift of 
electrons that scatter from the side dot on their path be- 
tween the two leads rcnormalizes the bare dot-2 phase 
shift 7722 ^ —tt/2 to produce an 7722 that jumps from —tt 
to at the p-h point. On moving away from £1 = — C/i/2, 
dot 1 gradually enters its mixed- valence regime, where 
there is no Kondo resonance and the T = conductance 
rises towards its unitary limit g = go- With increasing 
A, the Kondo state in dot 1 becomes more robust (as 
evidenced^^ by its larger Tk), pushing this upturn in g 
to larger values of \ei + Ui/2\. 

Raising the temperature progressively destroys the 
Kondo resonance and thereby increases the conductance. 
For fixed T > 0, 5 vs ei reaches a peak at ei — — [/i/2, 
where Tk is smallest and Kondo scattering is weakest. 
Figure [3ja) illustrates this behavior at three tempera- 
tures for £2 = 0, A = 0.03D. The double-dip structure 
surrounding the peak in g vs ei is qualitatively similar to 



the QPT feature in Fig. [2ljb). However, the temperature 
variation is very different in the two cases. In Fig. ^h), 
the decrease with increasing T of the conductance both 
at the peak and at the minima on either side is character- 
istic of the QPT. By contrast, conductance peaks arising 
for A(0) 7^ exhibit an increase with T of the extremal 
g values, as seen in Fig.[3Ka). 

To conclude, we have studied the linear conductance 
through a class of quantum-dot devices that can be de- 
scribed by a single Anderson impurity coupled to a con- 
duction band via a nonconstant hybridization function. 
Such devices can be tuned to a quantum phase transi- 
tion, marked by a near-unitary peak in the linear con- 
ductance that becomes more pronounced with increasing 
temperatures. The details of its evolution with temper- 
ature differentiate this conductance signature from simi- 
lar features arising from interference effects unrelated to 
quantum criticality. Our results demonstrate that these 
quantum-dot devices offer many advantages for the con- 
trolled experimental investigation of a rich array of many- 
body physics. 

We thank C. Lewenkopf, C. Biisser, and E. Vernek 
for valuable discussions, and support under NSF-DMR 
grants 0312939 and 0710540 (Florida), 0336431, 0304314 
and 0710581 (Ohio), and 0706020 (Tennessee). 



^ S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, U.K., 1999). 
^ M. Vojta, Phil. Mag. 86, 1807 (2006). 
^ A. Husmann et ai, Science 274, 1874 (1996). 
'' S. Sachdev, Rev. Mod. Phys. 75, 913 (2003). 
^ N. Mason and A. Kapitulnik, Phys. Rev. Lett. 82, 5341 

(1999). 
^ S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Sliahar, 

Rev. Mod. Phys. 69, 315 (1997). 
"^ D. Goldhaber-Gordon et al, Nature (London) 391, 156 

(1998). 
* H. Jeong, A. M. Chang, and M. R. Melloch, Science 293, 

2221 (2001); N. J. Craig et al, tbid. 304, 565 (2004); J. 

C. Chen, A. M. Chang, and M. R. Melloch, Phys. Rev. 

Lett. 92, 176801 (2004); A. Fuhrer, T. Ihn, K. Ensslin, W. 

Wegscheider, and M. Bichler, ihid. 93, 176803 (2004); R. 

Leturcq et a/., ihid. 95, 126603 (2005). 
® W. Hofstetter and H. Schoeller, Phys. Rev. Lett. 88, 

016803 (2001); M. Pustilnik, L. Borda, L. L Glazman, 

and J. von Delft, Phys. Rev. B 69, 115316 (2004); M. 

R. Galpin, D. E. Logan, and H. R. Krishnamurthy, Phys. 

Rev. Lett. 94, 186406 (2005); G. Zarand, C. H. Chung, P 

Simon, and M. Vojta, ihid. 97, 166802 (2006); R. Zitko 

and J. Bonca, Phys. Rev. B 76, 241305(R) (2007). 
"' R. M. Potok et al, Nature (London) 446, 167 (2007). 
" N. Roch et al, Nature 453, 633 (2008). 
^2 (a) L. G. G. V. Bias da Silva, N. P. Sandler, K. Ingersent, 

and S. E. Ulloa, Phys. Rev. Lett. 97, 096603 (2006); (b) 

ihid. 99, 209702 (2007). 



^^ The assumptions of a noninteracting dot 2 and a wide band 
simplify the analysis, but essentially the same properties 
arise from the full solution of Eq. ([T} for a weakly inter- 
acting dot 2, as will be discussed elsewhere. ^'^ 

" R. Bulla, T. A. Costi, and T. Pruschke, Rev. Mod. Phys. 
80, 395 (2008). 

^^ D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 
(1990); C. Gonzalez-Buxton and K. Ingersent, Phys. Rev. 
B 54, R15614 (1996); ihid. 57, 14254 (1998); R. Bulla, Th. 
Pruschke, and A. C. Hewson, J. Phys. Condens. Matter 9, 
10463 (1997); M. Vojta and R. Bulla, Phys. Rev. B 65, 
014511 (2001); L. Fritz and M. Vojta, ihid. 70, 214427 
(2004). 

^^ L. Vaugier, A. A. Aligia, and A. M. Lobos, Phys. Rev. 
Lett. 99, 209701 (2007). 

^■^ W. B. Lane, K. Ingersent, L. G. G. V. Bias da Silva, N. P. 
Sandler, and S. E. Ulloa, in preparation. 

^* K. Kang, S. Y. Cho, J.-J. Kim, and S.C. Shin, Phys. Rev. 
B 63, 113304 (2001); V. M. Apel et al, Microelectr. J. 
34, 729 (2003); C. A. Biisser et al, Phys. Rev. B 70, 
245303 (2004); P. S. Cornaglia and D. R. Grempel, ihid. 
71, 075305 (2005); P. Simon, J. Salomez, and D. Feinberg, 
ihid. 73, 205325 (2006); R. Zitko and J. Bonca, ibid. 73, 
035332 (2006); P. A. Orellana, G. A. Lara, and E. V. Anda, 
ihid. 74, 193315 (2006); A. C. Seridonio, M. Yoshida, and 
L. N. Oliveira, arXiv:cond-mat/0701529 (2007). 

1^ U. Gerland, J. von Delft, T. A. Costi, and Y. Greg, Phys. 
Rev. Lett. 84, 3710 (2000). 



